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ABSTRACT 

We present spectroscopic data for 180 red giant branch stars in the isolated dwarf irregular galaxy 
WLM. Observations of the Calcium II triplet lines in spectra of RGB stars covering the entire galaxy 
were obtained with FORS2 at the VLT and DEIMOS on Keck II allowing us to derive velocities, 
metallicities, and ages for the stars. With accompanying photometric and radio data we have measured 
the structural parameters of the stellar and gaseous populations over the full galaxy. The stellar 
populations show an intrinsically thick configuration with 0.39 < q$ < 0.57. The stellar rotation in 
WLM is measured to be 17 ± 1 km s , however the ratio of rotation to pressure support for the 
stars is V/cr ~ 1, in contrast to the gas whose ratio is seven times larger. This, along with the 
structural data and alignment of the kinematic and photometric axes, suggests we are viewing WLM 
as a highly inclined oblate spheroid. Stellar rotation curves, corrected for asymmetric drift, are used 
to compute a dynamical mass of 4.3 ± 0.3 x 10 8 M Q at the half light radius (r^ — 1656 ± 49 pc). The 
stellar velocity dispersion increases with stellar age in a manner consistent with giant molecular cloud 
and substructure interactions producing the heating in WLM. Coupled with WLM's isolation, this 
suggests that the extended vertical structure of its stellar and gaseous components and increase in 
stellar velocity dispersion with age are due to internal feedback, rather than tidally driven evolution. 
These represent some of the first observational results from an isolated Local Group dwarf galaxy which 
can offer important constraints on how strongly internal feedback and secular processes modulate SF 
and dynamical evolution in low mass isolated objects. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: dwarf — galaxies: individ- 
ual (WLM) — galaxies: kinematics and dynamics 



1. INTRODUCTION 

In models of hierarchical structure formation such as 
the current ACDM cosmologies, dwarf galaxies are fun- 
damental components of the mass assembly history of 
larger galaxies, such as our own Milky Way (MW). With 
dynamical masses below 10 10 M these objects offer im- 
portant observational tests for models, as they represent 



may be built (e.g., 


Navarro et al.|1997 


Moore et al.|1999 


Madau et al. 2001 


). Analyzing the survival of these low 



tant to constraining galaxy formation m odels ( Ricotti fc 
Gnedin|[2005| jGnedin fc Kravtsov||2006|) . Understanding' 
the factors that allow low mass dwarf galaxies to sur- 



vive to the present day may help explain discrepancies 
between observed and predicted distribution of subhalos 
around the MW. This requires an understanding of how 
internal and environmental effects are expected to shape 
the evolution of low mass systems like dwarf galaxies. 

With much shallower potential wells and lower metal- 
licities than their higher mass counterparts, dwarf galax- 
ies offer crucial laboratories in understanding how star 
formation (SF) proceeds throughout the lifetime of low 
mass objects. With total masses in the range 10 7 — 10 10 
M Q , dwarf galaxies have been shown to exhibit star for- 
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mation efficiencies much low er than higher mass galaxies 
( |Roychowdhury et al.l2009[ ). 

Environmental feedback offers one possibility for in- 
fluencing the stellar mass assembly rate. For exam- 
ple, is there a minimum halo mass which determines a 
galaxy's likeli hood of retaining baryons through reion- 



isation (e.g., Navarro et al. 1997| |Madau et al.||2001[ 
Ricotti fc Gnedin||ij()(]5| |Bov'll k Ricot,t,i||ij(J10| |Sawala| 



et al.||201ip 7 The apparent morphology-density relation 



in the Local Group, whereby more gas rich dwarf irregu- 
lars (dlrrs) are found at larger dista nces than the closer 



gas poor dwarf sphero idals (dSphs) ( Einasto et al.|[l974 
van den Bergh||1999|), is often invoked as evidence that 



environmental feedback has played an active role in low 
mass dwarfs. In this case both tidal and ram pressure 
stripping by the Milky Way could drastically alter the gas 
content and structure of dwarf galaxies, certainly playing 
a role in h ow their stellar content evolves over a cosmic 
time (e.g., Mayer et al. 2001b Kazantzidis et al. 2011| 



Lokas et al 



2010) 



Internal feedback such as SF and supernova (SN) 
driven winds may also play a role in modulating the stel- 
lar mass buildup of dwarf galaxies. In that case, due to 
the low potential well of the dwarf galaxies, the gas may 
not be retained during large episodes of s tar formation, 



short ening the SF lifetime in some cases (Dekel & Woo 



20031. This has been shown to have effects on not just 
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the chemical enri chment history, but also the structure of 
low mass dwarfs (Stinson et al.||2009| |Sotnikova k Rodi 



onov||2003 Governato et al.||2"Ul0p . With or without gas 
loss internal heating may regulat e the SF efficiency, lead- 
ing to lower star f ormation rates ( Brooks et al.|2007j ) and 
thicker structure (Kaufmann et ai.||2007p. Alternatively. 



low star formation efficiencies in dwarf galaxies may be 



due to the dependency of SF on H2. Recent work 1 


Bigiel 


et al.||2008 


Leroy et al. 


2008 


Schruba et al.||2011 


1 indi- 



lated with cold molecular hydrogen, H 2 , rather than the 
total gas surface density. Lower mass systems with low 
gas column densities do not allow self shielding of the 
H 2 , resulting in dissociation by the local radiatio n field 
( |Robertson k Kravtsov||2008| |Kuhlen et aL|2011| ). 



Further suppression of H2 formation is expected for 
low ISM metallicites, where fewer sites for dust forma- 
tion and longer cooling tim es for the gas would be present 
( Krumholz k Dekel|2011 ). If the SF efficiencies in dwarf 
galaxies are lowered due to their intrinsic low mass and 
metal content (both of which hinder the gas reaching the 
molecular phase), the amount of feedback need not be 
as large in order to e xplain the extended low level of SF 
( Kuhlen et al.| 2011 ). However the constraints on what 
relative contribution feedback still must play is not well 
known currently. Establishing whether low star forma- 
tion rates in dwarf galaxies are due to feedback or H2 
regulated mechanisms can be answered in part through 
analysis of nearby dwarf galaxies. 

Mechanisms that shape the evolution of the dwarf 
galaxy, whether internal or external, may impart obser- 
vational signatures onto the stellar and gaseous popula- 
tions in the galaxies. For example, the amount of pres- 
sure versus rotational support i s often tracked in simula - 
tions of dwarf galaxy evolution |k azantzi dis et al.|2011[). 



In the simu l ations of Ma yer et al 



et al. (2011 ); Lokas et al 



Hiiks (T<V/a < 5) are 



( |21X)lb| );|Ka 

rrs with thir 



Kazantzidis| 
in cold gas 



]2010|),d 

idally perturbed upon pericen 
a process invoked to explain the 



tre passage of the MW 
gas loss and change in stellar orbits that would be nec- 
essary to transform them into dSphs (V/a < 0.5). How- 
ever alternative theories exist in which the progeni tors of 
dSphs form in ho tter, thicker disks or spheroids ( |Kauf- 
mann et al.|2007 ) - this may in turn influence the SF ef- 
ficiency as the molecular gas for mation rate would be in- 
fluenced by the disk morphology ( Robertson k Kravtsov 
2008 1. And recent theor etical (IGovemato etal 
and observational work jl Sanc hez-Janssen et al i 
Roychowdhury et al.|2010 ), suggest that significant rear 
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rangement and evolution of baryonic structure can occur 
through secular evolution - whether via feedback (SF, 



SNe) , stellar migration, or global disk instabi lities (Stin- 
son et al. 2009 Sotnikova k Rodionov 2003). However, 
to connect the enrichment and kinematic history of the 
galaxy over most of the age of the Universe (the last 
~ 12 Gyrs), spectroscopic data (of resolved RGB stars) 
is required. In particular the relative amount of rota- 
tion and velocity dispersion of the stellar populations will 
evolve over the lifetime of the dwarf galaxy, and can be 
traced as a function of time. These spectroscopic obser- 
vations characterize the changing dynamics of t he tracer 



populations over the lifetime of the dwarf (e.g., van der 



Marel 2006 1 . This can be combined with deep photomet 



the galaxy ( Irwin k Hatzidimitriou| [l995 ) . Additionally 
stellar density profiles yield an estimate of the baryonic 
concentration, crucial for understanding the mass profiles 



van den Bosch 2001). Finally, coarse metallicity infor 
mation on a global scale can reveal chemical gradients 
across the galaxy. 

Therefore with structural analyses from photometric 
data, and the kinematics and ages from spectroscopic 
observations, it is possible to build a picture of how the 
mass distribution and dynamics of the galaxy evolved. 
While both internal and environmental feedback may 
have clear signatures on the dynamics of stellar and 
gaseous populations in low mass dwarfs, the question re- 
mains to what degree one effect dominates for a given 
dwarf galaxy. In isolated dwarf galaxies any changes are 
more likely due to internal effects - offering a window 
into how effective secular processes are at transforming 
a galaxy in the absence of strong tides. Any limits which 
can be placed on how feedback operates in these remote 
dwarf galaxies will offer important constraints in simula- 
tions of SF in low mass dwarfs - giving an observational 
framework that can constrain models of internal feedback 
mechanisms. 

Constraining the relative amount of internal feedback 
driven evolution requires studying galaxies that are likely 
to have been isolated over most of their history, un- 
like the nearby dSphs that have been influenced by the 
MW. While dSph distances (< 250kpc) put them well 
within the reach of modern 8m class spectrographs and 
imagers, they have most likely been tidally influenced 
by the MW, complicating interpretation. To date the 
large distances (> 500kpc) of isolated dlrrs have pre- 
vented substantive spectroscopic surveys, although these 
are the prime candidates to unravel the nature versus 
nurture question. Stellar structure and dynamics of the 
popu lations in the LMC h ave been studied in the past 
(c.f., van der Marel 2006), but due to its larger mass, 
close distance to the MW, and interaction with the SMC 
it is not ideal for studies of internal feedback driven evo- 
lution in a low mass dwarf galaxy. Ideally we wish to 
study a dwarf galaxy at large galactocentric distance, in 
which the evolved stars will likely have had fewer tidal 
distortions. With this we can understand for the first 
time, the dynamical state of the stellar component in an 
isolated dlrr. This can help constrain the initial condi- 
tions us ed for dwarf galaxy tidal transformation scenar- 
ios (i.e., Mayer et al.|2006"| ), which ma y be important for 
the re lative timescales of those models ( Kazantzidis et al. 
20TTI 



Galaxies such as the Wolf-Lundmark-M elotte dlrr 
(WLM, DDO 221; |Wolf|1910l|Melotte||1926[ ), sit at large 
galactocentric distances, and their complex stellar and 
gaseous populations trace how their structure and dy- 
namics evolve over a cosmic lifetime, with minimal exter- 
nal environmental influ ence. Observations of the bright 



supergiant population (iVenn et al. 
|2006| jUrbaneja et al.||2008t and 



20041 Brcsoli n et al. 
SErTSkillman et al. 



pgUl|Hodge k Miller l |'i995j [l^ee et al.|2005[) sampled the 
young populations of WLM, but ottered little insight into 
the earlier epochs of form ation and evolution. Similarly 
the photometric studies (Abies k Abies 1977 Minniti 
k Zijlstra||1997l |Hodge et al.||l<M)9l |Battinelli k bemers 



ric imaging studies, which trace the current structure of 



[jlstra||19971 |H 
2004||McConnachie et al.|2005||Jackson et al.|2007p were 
only able to provide global views of the evolved popula- 
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tion and were subject to degeneracies in age and metal- 
licity. We seek to answer two main questions on dwarf 
galaxy evolution: (1) what is the role of internal vs. ex- 
ternal feedback in shaping the structural and dynamical 
properties of the dwarf? and (2) what modulates the 
star formation efficiencies in low metallicity low mass 
systems - SNe feedback or H 2 regulated SF modes? To 
answer these, there needs to be a spectroscopic survey 
of RGB stars spanning 10 Gyrs in age, which produces 
abundances a nd velocities in an is olated dwarf galaxy 
like WLM. In |Leaman et al.|p009| (hereafter Paper I), 
we presented spectra tor 78 stars from FORS2 on VLT, 
which represented the first medium scale resolved spec- 
troscopic study of an isolated dwarf galaxy in the Local 
Group. In this paper we more than double our sample 
size, allowing for the first time a chance to study these 
questions in a truly isolated dwarf galaxy. 

1.1. Quantifying WLM's Isolation 

It is useful at this point to discuss the environment 
that WLM occupies within the Local Group, as the goal 
of this paper is to minimize the relative importance of ex- 
ternal (tides, ram pressure) effects in order to learn more 
about the possible impact of internal feedback effects. 

Derived and adopted properties for WLM are shown 
in Table 1 for clarity . WLM lies approximately 1 Mpc 
(iGorski et all 120111) from both the MW and M31, 



with its nearest neighbo ur being the sm all dSph Cetus 
(M ha if ~ 9 x 10 7 M^; |Walker et~aTJ |2009|), at 250 



kpc away ( W hiting et alT"1999), shown m Figure fl 



Karachentsev J 2005 ) calculate a tidal index for WLM o 
6 = 0.3, where they define 6 = max[\og(M k /Df k )} + C 
to be the amount a galaxy is acted on by its largest 
tidal disturber^] . With the exception of Tucana (-0.2), 
Pegasus (-0.1), Aquarius (-0.1), and Leo A (0.2), WLM 
is one of the five least tidally disturbed and most isolated 
galaxies within a ~ 1 Mpc sphere of the MW. Projection 
of WLM's heliocentric velocity (v sys ~ —130 km s -1 ) 
towards the Local Group centre of mass results in a 
velocity with respect to the barycentre of vlq = —32 km 
s _1 , implying that it has just passed apocentre. With 
its current distance and that velocity, and an inferred 
Local Group mass of 5.6 x 1O 12 M (based on M31-MW 
orbit timing argumen ts which assume an age of 13.7 Gyr 
for the universe; c.f., Lynden-Bell 1981), the maximum 
apocentre that WLM could have is approximately 1.3 
Mpc. Assuming a completely radial orbit, the implied 
orbital period for WLM to reach its pericentre with the 
Milky Way is 11 - 17 Gyr. Thus WLM, in addition 
to being currently quite isolated, has had at most one 
pericentre passage in its lifetime (which would have 
been at least 11 Gyr ago at z ~ 2.5), meaning its total 
evolution has been much less dominated by tidal and 
ram pressure effects from the Milky Way than other 
Local Group galaxies. 

The paper will be presented as follows: in §2 we 
will present the data sets with a focus on the obser- 
vations and reductions. §3 will discuss the spectral 

1 As noted in that work, is the fcth disturber galaxy at a 
distance D^, and © = corresponds to an object with a Keplerian 
orbital period about MW/M31 equal to 1/H. For reference, Sgr 
dSph has © = 5.6 



analysis techniques used to reduce and analyze the data 
in a homogeneous way to extract the chemodynamical 
properties for this study. §4 will present results on 
the structure and kinematics of the stellar and gaseous 
populations, §5 a discussion of the mass estimates, and 
§6 will focus on the evolution of the stellar dynamics. 

2. OBSERVATIONS AND DATA REDUCTIONS 

The data presented in this paper comes from observa- 
tional campaigns on several instruments an d telescopes. 
78 stars w ere observed with the FORS2 (Appenzeller 
et al.|1998 ) spectrograph on V LT and this data was pre- 
sented, in Leaman et al. ( |2009 ). These spectroscopic ob- 
servations are supplemented here with spectra of 140 new 
stars observed at h igher resolution with the DEIMOS 
(Faber et al. 2003) spectrograph on Keck II . V and 
i band photometry from the INT Wide Field Camera 
covering a 36 x 36 field of view was used for several 
analysis steps. This pho tometric data is presented in 



McConnachie et al. ( 2005 1 and has been converted to V 
and 1 using the calibration in that work. Near infrared 
JHK data from the UKIR T/WFCAM telescope was also 
used (Tatton et al.||2010). I n addition we use VLA radio 
observations ol WLM from Kepley et al. (2007) which 
the authors kindly made available for use here (A. Ke- 
pley, priv. comm.). Figure illustrates the coverage of 
our slitmasks for both the FORS2 and DEIMOS obser- 
vations. 

The new DEIMOS observations were taken during 3 
nights of September of 2009. Two 16 x 5 slitmasks 
covering the body of WLM were observed (see Fig. [2| 
during seeing that ranged between 0.7 — 1.2 , with a 
median value of 0.8 . The instrument setup used the 
12001/mm grating with OG550 blocking filter, which 
yielded a spectral resolution of ~ 1.4A through 1.0 
slits. The spectral range spanned roughly 7800 — 9300A. 
Wavelength calibration was aided with the use of NeAr- 
KrXe arc lamp exposures, in addition to the standard, 
quartz flat field calibrating exposures. Each mask had 
10 x 30 minute exposures plus 3 x 20 minute exposures, 
all taken during good seeing and weather conditions, 
yielding 6 hours integration time per slit. 

As the DEEP2 reduction pipeline is not optimized for 
faint stellar spectra, the data were pipeline processed us- 
ing a privately written reduction procedure (R. Ibata). 
In brief the pipeline performs the requisite calibration 
tasks to characterize the slit position on the CCD and 
remove any instrumental signatures, and non-orthogonal 
projection biases. Details of the reductio n procedure 
used f or the DEIMOS spectra is described in [Ibata et al.| 
(20111 and references therein. Extraction of the spec- 



tra m each slit results in a single stellar source which 
is further subject to wavelength calibration and contin- 
uum normalized, before being radial velocity corrected 
via cross correlation against a suite of template spec- 
tra. Due to the extreme faintness of the stars, even 
10m class telescopes require coaddition of 6-8 hours of 
integration time to produce S/N > lOA^ 1 data appro- 
priate for resolved chemodynamic studies such as this. 
Due to concerns about the poor performance of the flex- 
ure compensation system on DEIMOS at the time of 
our observations, the spectra were extracted individually 



4 



Leaman et al. 



and processed as single data products and coadded later. 
The signal to noise per angstrom for the FORS2 spectra 
ranged from 17 < S/N < 30, and 6 < S/N < 30 for 
the DEIMOS data. Figure [3] shows three representative 
spectra of stars observed with our DEIMOS configura- 
tion. 

In addition to the two masks on WLM three galactic 
clusters were observed for calibration purposes - the old 
open cluster NGC 6791, and globular clusters Pal 14 and 
NGC 7078. Single mask setups yielded 5 — 50 stars per 
cluster, which were observed 2-3 times per night through- 
out the observing run. This allowed us to check the pre- 
cision of our metallicities and assure that we are bringing 
the [Fe/H] estimates onto a common, calibrated scale. 

3. SPECTRAL ANALYSIS 

3.1. Radial Velocity Measurements 

Radial velocities were measured from the three strong 
calcium II triplet (CaT) lines (A ~ 8498, 8542, 8662A) 
which in the case of the DEIMOS data, had a veloc- 
ity correction applied from the sky OH lines in order to 
mitigate errors due to the instrument rotation. The av- 
erage rms of the wavelength solution for the blue and 
red side of the DEIMOS chips were 0.25 and 0.27 km 
s _1 respectively. For the FORS2 spectra, the low sig- 
nal to noise of the individual frames necessitated that 
we perform cross correlation radial velocity calculations 
on the combined spectra, rather than on each individ- 
ual image. As such heliocentric velocity corrections were 
tailored to the individual exposures and applied prior to 
combining the spectra, due to the long temporal base- 
line (roughly four months) of the FORS2 observations. 
Once shifted and combined, the spectra were ready for 
radial velocity computation with the aid of template stars 
and a Fourier cross correlation routine (fxcor) . For the 
FORS2 data, a total of 23 template radial velocity stars 
observed with the same instrument setup were used with 
the cross correlation routine. In the case of the DEIMOS 
data, the stellar spectra were cross correlated against a 
single synthetic template around the CaT region. Both 
computations provided in-line error estimates, with the 
velocity errors ranging from 1.0 < SVhei < 8.0 km s _1 
with a mean of (SV he i) = ±2.3 km s" 1 for the DEIMOS 
stars, and a range of 3.0 < SVhei < 10.0 km s _1 with 
a mean of (SVhei) — ±5.0 km s _1 for the FORS2 stars. 
For the FORS2 reductions systematic velocity errors due 
to a star's position in the slit were removed by centroid- 
ing the stars relative to the slit centre. The fact that 
this procedure was done on combined spectra resulted 
in small absolute corrections, as the i/n statistics meant 
that the individual slit errors were minimized in the com- 
bination and correction steps. The typical shift for the 
slit error on an individual exposure is approximately 6- 
9 km s _1 for the FORS2 stars, and we note that this 
shift produces negligible uncertainties in the equivalent 
width error. The final average absolute corrections to 
the slit centering errors on the combined spectra were on 
the order of < 1.5 km s" 1 for the FORS2 stars. 

3.2. Membership Determination 
Due to its distanc e off the plane of th e Milky Way 



WLM. Nevertheless, the DEIMOS sample of WLM cov- 
ers a much larger field of view than the FORS2 sample 
presented in Paper I, which means the chance of fore- 
ground stars in the outer regions is increased. Addition- 
ally, as we are possibly sampling member stars at larger 
radii, they stand to be at larger velocities with respect to 
the central regions. Removing contaminants in a galaxy 
that is known to be rotating is not trivial, and in this case 
we culled the sample with an iterative maximum likeli- 
hood method. Rather than a global 3— sigma clipping, 
we adopted a spatially binned clipping routine which was 
much more robust, and by binning along the major axis 
avoided artificially inflating the o~ v estimates. This proce- 
d ure was based on the .j oint maxi mum likelihood method 
oflWalker et al.| |2006j ) (see also |Gunn fc Grrffin]p79" 



Hargrcaves et al. 1994 ) , which we extended here to spa- 

^ — i — r-> — ' — i ' h p i " ■ i i i i i i 



tial bins of 25 stars each, so as to accurately throw out 
contaminants when the rotational profile enhances the 
spread of velocities at different positions in the galaxy. 
The algorithm simultaneously determines the dispersion 
(<j v ) and average velocity ((it)) in a given bin by maxi- 
mizing the probability: 



1 N 

i=l 



1 N I 



H-l+^y: VV.^' -y Mar). 

(i) 

With this we culled the sample to 180 total member stars, 
with only one of the previous FORS2 stars removed. Fig- 
ure [4] shows the distribution of probable member and 
contaminant stars in our catalogue. In Figure [5] we plot 
the 2-D projected map of WLM with the HI gas con- 
tours and stellar sources determined to be members in 
our sample, both colour coded by velocity. 

3.3. Spectroscopic Metallicities 

As in Paper I, [Fc/H] derivations were done through 
the empirically calibrated Calc ium II triplet method 
(Armandroff & Da Costa 1991) whereby the summed 
equivalent width of these lines correlates with [Fe /HI on 



a scal e calibrated to galactic globular clusters ( Cole et al. 
2004|. We here use the recently updated calibration 



given in Starkenburg et al. (2010). The details of the 
equivalent width measurements and placement onto 
the metallicity scale will be discussed in a forthcoming 
partner paper (Leaman et al., 2011b) which will deal 
with the chemical analysis and interpretation of WLM. 
Of our 180 radial velocity members we have derived 
[Fc/H] measurements for 126 of the stars which had high 
enough signal to noise (> 10A _1 ). The uncertainties 
were estimated in the fashion described in Paper I, 
with the average internal error on metallicity being 



A 



Fc/H] 



±0.26 dex. 



(I, b = 75.85, -73.63; |Gallouet et"aT|l975| ), the contami- 
nant fraction is expected to be very low m the direction of 



3.3.1. Repeat Observations 

While estimates of our errors exist for both [Fe/H] and 
velocity measurements, more information about how pre- 
cise these measurements are can be gleaned from obser- 
vations of stars that are in common between the FORS2 
and DEIMOS dataset. We ended up with 6 stars that 
overlapped between the sample for which to check ve- 
locities - and only 2 stars for which repeated [Fe/H] 
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measurements can be made. The values for these re- 
peat measurements are shown in Table 2 . The mean 
offset ((X) = jjH(xi — Xj)) and scaled median absolute 
deviation of the offsets (S = 1.48 x med(\xi — med(xj)\)) 
are X\p e /m — 0.24; S\p e /m = 0.02 dex, and Xy e i = 
— 2.5; Svei = 13.9 km s _1 respectively. Given the ve- 
locity and metallicity uncertainties discussed in the pre- 
ceding section, these repeat observations are consistent 
within the errors. 

3.4. Age Derivations 

At the distance of WLM, deriving ages for individual 
stars is not possible with precision greater than ~ 50%, 
and in reality is simply providing a coarse correction to 
the assumption that metallicity and colour are perfectly 
correlated. The reader is referred to Paper I for the dis- 
cussion of the age derivations in detail, but relevant up- 
dates are briefly outlined here. The well calibrated ho- 



mogeneous photometry from McConnachie et al. ( 2005 1 
is crucial in deriving accurat e ages, but th is method still 

Tolstoy|2003|) 

3ct 



suffers from difficulties (c.f. 



namely the 



possibility of internal differential reddening, and poorly 
constrained [a/Fe] values for the RGB stars as di scussed 
in Paper I. Given the extended SFH of WLM ( |Mateo 
1998 Dolphin|2000[ ) and the gas rich nature of dlrr galax- 
les, we expect a significant range in red giant ages in the 
dataset. 

With our current spectroscopic sample, [Fe/H] val- 
ues are now found < —2.5 dex, which necessitated us- 



ing th e Yale-Yonsei evolutionary tracks ( Demarque et al. 
2004 1 , which sample much lower metallicities. In Paper 
1 we provided age estimates based on both [a/Fe]= +0.3 
and solar. Here we use a compilation of high disper- 
sion spectroscopy of various stellar populations in WLM 
to constrain the mean [a/Fe] as a function of metal- 



licity. Colucci et al. (20101 observed the lone globular 



cluster (GC) of WLM7~WLM-1 to have [Ca/Fe] = 0.25 
and [Fe/H] = -1.71, which we adopt as the [a/Fe] values 
for all our targets with similar or lower metallicities. Of 
course GCs may undergo a separate formation and chem- 
ical evolution than field stars, however in the MW, LMC 
and Fornax the m etal poor clusters show similar [a/Fe] 



as the field stars (Pritzl et al. 2005 Hill 



1997 Letarte 



et al. 2006) , and Oa is not affec ted by GC mixing pro 



cesses (e.g., Gratton et al.||2004[) . Stars more metal rich 
than this are assigned an [a/FeJ based on a third order 
polynomial, that is fit to the WLM-1 GC value, as well 



as t he B, and A supergian t studies of Venn et al. ( 2003 1 



and |Bresolin et al.| (120061). The position ot each star on 
the colour magnitude diagram was interpolated among 
a grid of Yale-Yonsei isochrones at appropriate values of 
[Fe/H] and [a/Fe] for the star. 

The random errors on ages for individual stars are 
±50%. Systematic errors are present in these estimates 
due to the variable a-enhancement, extinction, and AGB 
contamination. We have estimated the magnitude of 
these systematic uncertainties in our ages using artifi- 
cial star tests, which will be described in more detail in 
the forthcoming paper (Leaman et al., 2011b). In the 
relevant figures in this work we have overlaid both the 
individual random error in a star's age, as well as the 
combined systematic uncertainties. 

4. RESULTS 



In this work we have added DEIMOS spectroscopic 
observations of 103 member stars to supplement the 77 
member stars from Paper I. We have determined the ra- 
dial velocity values for these 180 individual member RGB 
stars in WLM, of which 126 also have [Fe/H] and ages 
derived. These new objects allow us to sample out to a 
projected radius of 3.4 kpc (<~ 3.5 disk scale lengths) in 
WLM, assuming the adopted distance of 985 kpc. This 
spatially extended sample can be used to examine the 
structure, kinematics, and chemical evolution history of 
this galaxy for the first time in conjunction with derived 
age estimates. We now present an updated discussion of 
the structural and kinematic parameters. 

4.1. Structural Parameters 

To derive the structural parameters of WLM we used 
the photo metric V,I, and JHK catalogues, as well as the 
HI data of Kepley et al. (2007]). The colour magnitude 
diagrams are shown in Figure |6j including the INT V, 
I photometry as well as the infrared UKIRT / WFC AM 
photometry of WLM. We have isolated several popula- 
tion tracers, as shown by the coloured areas in the figures. 
For the V, I photometry a strict RGB locus is shown by 
the red box. Th e tip of the RGB loca t ion wa s taken from 
the analysis of McConnachie et al. (2005) which used 
the same photometric catalogue as our present work. 
This selection represents the least contaminated red giant 
branch population in WLM. Outside of this lie primarily 
excessively young or metal poor RGB and luminous AGB 
stars. In the JHK diagram the oxygen rich M-stars (light 
blue) and the young (~ 1 Gyr) carbon rich C-stars are 
identified, as those AGB stars above the TRGB and to 
the appropriate side of the metallic ity dependent colou r 
cut found for WLM in the work of jTatton et al | ( |2010 ). 
The stellar density profiles for the full V,l (black), strict 
RGB (red), and C-star (blue) populations are shown in 
Figure J7J constructed by binning stars in geometrical 
radii which have been referred to an average ellipticity of 
e = 0.55 at a position angle of 179 degrees. The inner re- 
gions of the two INT populations suffer from crowding, 
and so the subsequent fits of the surface density pro- 
files exclude these regions. Four p rofiles fi ts have been 
applied; Plummer (Plummer 1911), King ( |King| [l962), 
de Vaucouleurs and exponential profiles. For the rest of 
this paper we adopt the best fitting exponential surface 
brightness profiles, £(i?) = T,Qe~ R ^ Rd , and convert the 
exponential scale length, Rd — 3.45 (for the full V,I sam- 
ple), into the 2D projected half light radius, r^ = 5.78 , 
and where necessary the 3D depr qjected half light r adius, 



Wolfetal. 



(2010). The 



r/j 3 = 7.68 , using the relations of 1 
background level in each population was simultaneously 
fit as a constant, and was found to be consistent with 
selecting out by eye isolated regions of our field of view 
to estimate the background level. 

The 2D projected spatial distributions of these three 
populations and the HI surface density are shown in Fig- 
ure [HI The mean and standard deviation of the back- 
ground level was determined from regions at the outer 
limits of our field of view. Contour levels start at 2a 
above the background and each subsequent contour level 
is n = 1.2 times higher than the last. The HI radial 
extent is not signifi cantly larger than th e stellar compo- 



nent as reported in |Kepley et al.| ([2007]) and the recent 



G 
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single dish observations of WLM by Hunter et al. (2011). 
The latter work also noted that the density profile of the 
HI shows a smooth drop off to the detection limit, with 
no sharp truncation. The strict RGB population shows 
a slightly more spherical distribution than the HI, but 
otherwise the general shape of the evolved stellar popula- 
tions are consistent with the HI population. As expected 
the young C-star population occupies the central regions 
of the gal axy, where the most recent star formation has 
occurred (Dolphin 2000). These centrally concentrated 
younger populations are similar to what has been found 
in other Local Group dwarf galaxies (e.g., [Battaglia et al. 
2006| |Hidalgo et al.|2009||de Boer et al |201l[ ). However 
we note that with characteristic ages of ~ 1 Gyr, these 
stars would still have experienced ~ 3 rotation periods. 
Ellipticity and photometric position angle, (PA p hot), are 
calculated as a function of r adius by fitting ell ipses us- 
ing the task MPFITELLIPSE ( |Markwardt]j2009[ ), on the 
contours at various radii for each population. 

The resulting position angles and ellipticities are 
shown as a function of radius in Figure [9] There is 
excellent agreement between the HI and stellar struc- 
tural properties with only the evolved RGB population 
showing a slightly rounder ellipticity. This is consistent 
with literature values of past photometric stu dies of 



(Abies k Abies 1977 Minniti k Zijlstra 19971. The 



close agreement m position angle and ellipticity of the 
stellar populations in WLM is in contrast to some lower 
angular momentum dSphs such as Phoenix and Forna x 



( Martmez-Delgado et al. 
WLJV1 clearly is flat tenec 



19991 IBattaglia et al. 2006 ) 



more than the nearby dSphs 



in the Local Group (Mateo|1998l - which may be due to 



the strong rotation as we shall see in the next section 



4.2. Internal Kinematics 

With nearly equal rotation and dispersion in the stellar 
component, care must be taken in interpreting WLM's 
velocity data. In rotating systems, classical mass esti- 
mates rely on the HI gas, which has negligible random 
motion, making for relatively simple analyses equating 
the enclosed mass to some circular velocity. The under- 
lying mass profile can also be ascertained through veloc- 
ity dispersion profile b ased methods - such a s spherical 
Jeans modelling (e.g., [Battagli a et al.| [2008 ) . In each 
of these cases, only one type of dominant kinematics is 
present in the tracer population - rotation or dispersion. 
The challenge in our case is to convert WLM's composite 
dynamical configuration into a singular rotating or dis- 
persion dominated dataset. Without these corrections 
a number of systematic biases arise when interpreting 
rotational and pressure supported velocity distributions 
- such as encountered whe n analyzing an obj ect in the 
configuration of WLM (c.f. Lokas et "aL]|2010 |. Detailed 
dynamical modelling will be presented m a forthcoming 
paper (Leaman k Battaglia, in prep.). However the anal- 
ysis of rotation and velocity dispersion profiles used in 
the remainder of this work requires accurate removal of 
the biases that mixed rotation/dispersion systems show, 
which we discuss now. 



4.2.1. Stellar Velocity Field 

As discussed in Paper I, the stellar velocities are 
markedly different in WLM from most of the dSph galax- 



ies. The stars in WLM of all ages/metallicities appear 
to be rotating, in contrast to the purely random stellar 
motions in Local Group dSphs. We begin by looking at 
the velocity structure of the stellar and gaseous popu- 
lations of WLM. As shown in Figure [TUJ WLM's stellar 
velocities, while showing rotation greater than seen in the 
dSphs of the Local Group, rotate with half the velocity 
of the HI gas, and lag behind the HI by ~ 15 km s _1 . 
Clearly the stellar rotation is decoupled from that of the 
HI, with the suppressed rotation velocity of the stars due 
to the effect of asymmetric drift as we will demonstrate 
later. In Figure [TTJ we show the projection of the rota- 
tion velocity of the stars as a function of geometric radii 
(again referred to an ellipticity of e = 0.55). The binned 
data points are fit with an isothermal sphere rotation 
curve of the form V (r) 2 = ATra$a\(l — ai/r)(atan(r/ai)) 
( Kepley et al.|200"7 ), which slowly rises to a value of ~ 17 
km s _i at the last measured point. 

As later we wish to derive a line of sight velocity dis- 
persion profile, it is important that we accurately charac- 
terize the velocity field of WLM. We therefore explored 
different models to characterize the rotation of the stars 
in WLM. A few methods that were tested were: fitting 
the stellar velocity field simply along a strip of the major 
axis, describing the rotation field by iterating through 
fits at differing position angles, or treating the stars with 
the classical "spider diagram" velocity field: 

v(x,y) = v sys +v rot {R)cos(9)sin(i)+v rad (R)sin(9)sin(i), 

... ( 2 ) 

where 9 is the azimuthal angle and i the inclination of 
the galaxy. 

The most consistent determination, as judged by the 
smallest deviation at all radii in the rotation-subtracted 
velocities, was a simple linear fit along the velocity pat- 
tern of the major axis. The rotation subtracted veloci- 
ties are shown in the second panel of Figure |11[ with the 
mean values of the stars staying close to the systemic 
velocity at all radii. The velocity of the stars off the ma- 
jor axis shows little deviation from the model, suggesting 
that the stars are rotating cylindrically - however due to 
the uncertainties in the inclination and disk thickness of 
WLM we can not examine this in detail. 

4.2.2. Velocity Dispersion Profile 

Deriving a line of sight velocity dispersion (a v ) profile 
for these stars requires removal of the rotation signa- 
ture, so that the a v profile is not artificially enhanced 
at large radii. As we showed above, WLM exhibits 
significant rotation, so correcting for this magnitude of 
rotation in a galaxy with this structure is not trivial. 
Using the best fit velocity field model described in the 
previous paragraph, the rotation was removed and we 
could proceed to construct a velocity dispersion profile 
that was unbiased at large radii. The dispersion was 
calculated applying the likelihood method of Equation 
1, using the rotation subtracted velocities. Errors on 
the dispersion estimates were taken to be the value at 
the 1 — a confidence interval of the likelihood parameter 
space. The resulting dispersion profile is shown in the 
third panel of Figure [TTJ as solid black dots. The open 
blue circles in that figure show the resulting a v profile 
if we were to compute the dispersion profile on the raw 
velocities that hadn't been rotation subtracted. The 
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artificial enhancement of the dispersion values at large 
radii is clearly evident due to the strong streaming 
motion of the stellar velocity field. The radial velocity 
dispersion profile for WLM does n ot appear to be flat , 
unlike what is seen in many dSphs ( | Walker et aLp009l ), 
however interestingly the dE NGC 147 which is of 
similar mass, luminosity, and shape as WLM shows a 
drop off in a v at larger radii as well (Geha et al. [2010). 
This and other nuances of the dispersion profile will 
be discussed in depth in a forthcoming paper on the 
detailed mass profile of WLM (Leaman & Battaglia, in 
prep.) 



11 



we show the 



In the bottom panel of Figure 
ratio of rotational to dispersion support m WLM as a 
function of radius, V/a. The dashed line represents the 
expected value of V/a for an oblate isotropic rotator at 
the eccentricity of WLM. The flattening of WLM then 
is most likely due to rotation, rather than anisotropy 
in the outer regions. The value of V/a ~ 1 within 
one effective radius is much less than the HI, which 
shows V/a ~ 7. Nearly equal pressure and rotation 
support for the stellar component is also much lower 
than what simulations typically use for progenito r dlrrs 
(i.e., |Mayer et al.|[2001a| Kazantzidis et al.| 2011[). The 



measured V/a may reflect some minor interactions 
between WLM and other dwarfs in the Local Group, 
however as discussed in the next section, the alignment 
of the kinematic and photometric axes of WLM rules 
out significant interactions. Therefore if this value of 
V/a = 1 for the isolated WLM dlrr is representative of 
typical dlrr progenitors, it would have implications for 
the efficiencies (and thus timescales) of dwarf galaxy 
tidal transformation scenarios in such simulations. 



4.3. Photometric and Kinematic Axes 

Information on the orientation and potential of WLM 
comes from comparing the position angle of the kine- 
matic (PAkm) and photometric (PA p hot) data. The kine- 
matic position angle defines the major axis of rotation of 
a galaxy, and in cases where recent mergers have oc- 
curred or the observer is viewing a triaxial system, may 
not always lie along the same axis as the photometric 
posi tion angle. To comp ute the PAkm we proceed as in 



e.g., Walker et al. (2006). For every star, or interpolated 



HI velocity point in the tangent (£, r)) plane, an angle of 
bisection is calculated based on the data point's position 
off the above determined photometric minor axis. The 
galaxy is split along this bisection line and velocity dif- 
ference is computed by calculating the average velocities 
above and below the line and subtracting them. The bi- 
section line is then rotated by one degree for the HI or 
to the next star and the proced ure is repeated, building 
up the curves shown in Figure |12| A hyperbolic tan- 
gent function was fit to the data and the position angle 
determined based on the bisection angle where the ve- 
locity difference is minimized. With a position angle of 
zero northwards in the tangent plane on the sky, the 
kinematic position angle is defined as 180 degrees + or - 
the angle of bisection where the minimum velocity differ- 
ence occurs. The uncertainty on the angle is determined 
by calculating the difference from the bisection line that 
would be produced if the velocity difference was changed 
by the average uncertainty in a stars velocity. We find 
consistent values for the kinematic position angle of the 



stellar and gaseous components of WLM, with the lat- 
ter values in excellen t ag reement with the HI studies of 



Jackson et al. (20041 and Kepley et al. (2007), who used 
other methodology. 

With the photometric and kinematic position angles 
calculated, it is instructive to look at their relative 
alignment, typically designated as the parameter ^ = 
\PAkin — PA p hot\- The quantity of ^ is sometimes used 
as a signature of mergers, where the structural and dy- 
namical misalignment would be high, or to give insight 
into non-circular gaseou s motions due to infall or outflow 
( |Fraternali et al.||2004[ ). For galaxies in the ATLAS 3D 
study ~ 90% of galaxies wer e found to have alignments 
< 15° ( Krajnovic et al.|2011 ). For the stellar component 
of WLM we derive ^ = 1 ± 3°, and W = 2 ± 2° for the 
HI. The fact that the alignment of the kinematic and 
photometric position angles is close, also constrains the 
viewing angl e and shape. As discussed recently by |Lokas| 
et al. (2010), when rotation is significant in a triaxial 



system viewed perpendicular to the longest dimension, 
rotation other than the major axis is usually evident, as 
the likelihood of the rotation being aligned with either 
axis is small in a potential that permits more irregular 
orbits. In contrast an axisymmetric system will show no 
rotation along the minor or intermediary position angles, 
thus the kinematic and photometric position angles will 
show strong alignment like we observe in WLM. Addi- 
tional constraints against WLM being triaxial come from 
its position on the V/a vs. e anisotropy diagram. As we 
shall see in the next section WLM is flattened mostly by 
rotation rather than anisotropy, consistent with the in- 
terpretation that we are viewing a highly inclined oblate 
spheroid. 

4.4. How Thick is WLM? 

From Figure [9] the ellipticity and position angle of 
the gas and stars allow us to characterize the projected 
distribution of those populations. To untangle the in- 
trinsic shape of WLM is more complicated. If we are 
observing an oblate spheroid that is rotationally flat- 
tened, it will have an intrinsic thickness, qo (for this 
orientation, q = c/a, with c being the short axis and 
a being the long axis) - which is difficult to disentan- 
gle without a constraint on the viewing inclination, i. 
The relation between the intrinsic and projected axial 
ratio is shown in Figure [13] for the range of inclinations 
de termined in the tilte clring analysis of the HI data 
in Kepley et al. ( |2007[ ). We note that because their 
inclinations were derived assuming a thin disk in the 
tilted ring analysis, they therefore provide lower limits 
on the possible inclination. Thus given the high incli- 
nation for WLM, the range of measured projected axial 
ratios (b/a) for the stellar and gaseous components (hor- 
izontal lines), suggest that the intrinsic thickness ranges 
between 0.38 < qo < 0.57. With the data presented 
here, we have evidence that WLM's vertical structure is 
quite extended in both its stellar and gaseous compo- 
ne nts. This is confirmed using the techniques described 



van der Marel et al.| ( |2002[ ), whereby WLM can be 



approximated as a flattened spheroid with axial ratio 



q{R) = {2- 



\ 2R E 
t R 



vWg) \2n-i/2 



a(R) 



) 2 ]) 1 ' 2 , assuming isotropic ve- 



locities for the stars (which are assumed to lie in a larger 
isothermal halo). Using the derived exponential scale 
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length from §4.1, the intrinsic axial ratio (q) for WLM 
calculated using this method again lies between 0.38 and 
0.56 at all radii - in good agreement with Figure [13] 

This value is consisten t with the recent study of |Roy-| 
chowdhury et al. ( |2010 ), who showed that the average 
intrinsic gas disk thickness of dwarf irregulars is a much 
higher (qo) — 0.6, instead of the 0.2 that is canonically 
assumed for stellar and gaseous disks in higher mass 
galaxies. Our values for the stellar thickness are also 
in line with studies of distant d warf galaxies that seem 
to show thicker stellar structure ( | Sanchez- J anssen et al. 



2010). 



We might ask how does this thick configuration 
compare to other galaxies where direct measurements of 
the vertical height are possible? While we are unable 
to determine whether the extended vertical structure 
in WLM is due to just a single component flattened 
spheroid, or multiple disks of varying thickness, we can 
compare the vertical extent to studies of edge on galaxies 
that exhibit thick disks. In Figure [l4j we plot the radial 
and vertical scale heights for W LMalong with a sample 
of edg e on spiral galaxies from |Yoachim fc Dalcanton| 
( 2006 [) , and the two large spirals in the Local Group (c.f. 
Collins et al. 2011). The values for WLM were derived 
by fitting the major axis stellar surface density with an 
exponential profile, and the minor axis with a profile 
of the form T,(z) oc sech 2 (z/z ), for which we derive a 
vertical scale height of zq — 705 ± 28 pc. Assuming 
an inclination of i = 79 we show the corrected vertical 
scale height marked on the figu re as the open circle. 
The b est fit lines are taken from lYoachim fc Dalcantonl 
(2006) and indicate that WLM has a similar relative 
vertical scale height as expected for its radial scale 
length. Direct confirmation of a distinct thick disk in 
WLM would only be possible after identifying a velocity 
lag or age difference in a statistical sense, which is 
beyond the scope of our spectroscopic sample. However 
the vertical structure of WLM is clearly extended to 
an amount consistent with other galaxies that show 
increased vertical scale heights. 

5. MASS ESTIMATES 

Examination of the velocities shows that WLM is 
clearly in a different structural and dynamical state than 
dSphs, however does that mean their mass estimates also 
differ? Detailed mass modelling can provide answers to 
both what the stellar velocities tell us about the past his- 
tory of WLM, and about the underlying potential. As 
we have both stellar and HI velocities there exists the 
opportunity for much more detailed modelling than is 
available from dSph data sets. With this we can com- 
pare how various mass modelling techniques perform on 
a system where there is strong rotation and flattening. 
Various dynamical models to describe the mass profile of 
WLM will be presented in a forthcoming paper but for 
now we can compute simple mass estimates in two ways. 

5.1. Rotationally Derived Mass 

We can use the inherent rotation of the stars to com- 
pute an estimate of the circular velocity of WLM, pro- 
vided we remove the component of random motions that 
are suppressing the rotation curve. As noted earlier, the 
stellar rotation curve differs in magnitude to the HI due 



to the effect of asymmetric drift. To directly compare 
the rotation curves of the stars and gas, and allow for in- 
dependent circular velocity estimates we must correct for 
this effect. To analyze the stellar kinematics in a purely 
rotational sense corrections due to asymmetric drift were 
applied. 



dp 
dR 



d In a\ 
dlnR 



(l--f-) + ifcr 



Here 
velocity, 



dz 
(3) 



is the inclination corrected observed rotation 
p is the stellar density, a the velocity disper- 
sion in a given component, and v c is the drift corrected 
rotation velocity. We foll ow the outline and assump- 
(20011: namely that the last term 

= 0.5, the density profile is of 

R 

-R/Rd), and the slope of the ra- 



tions of Hinz et al. 



can be neglected, - 

the form p cx exp( 
dial component of the velocity dispersion is flat in the 
outer regions of the disk. The equation then reduces to 
v c = v l + cl(2fi/-Rd - 1). We note that the assumed 
anisotropy is chosen primarily because it reproduces flat 
rotation curves in large galaxies, however we have little in 
the way of constraints on it at this time. We use a nom- 
inal inclination of i = 79° for the inclination corrected 
values of dispersion and rotation, allowing us to com- 
pute the corrected stellar circular velocity, v c . As noted 



earlier, the tilted ring analysis of Kepley et al. 
assumed a thin disk, therefore their inclination w 



1)2007 1 
tiich we 



adopt here should provide a conservative lo wer limit. The 
results of this correction are seen in Figure |10| The final 
inclination and asymmetric drift corrected curve shown 
by the blue binned points matches well with the HI gas 
- giving us confidence that our dispersion measurements 
are correct for the stars. 

With this corrected stellar rotation curve we have es- 
timates of the circular velocity at various projected radii 
in WLM, with which we can compute the associated ro- 

tationally derived mass as M(< r) = ^f. In Fi gure 15 
we show these dynamical mass estimates as a function or 
radius as the black squares. The black solid and dashed 
lines show the enclosed mass curves based on the HI ro- 



tation analysis from the study of Kepley et al. (2007), 
assuming an isothermal sphere model. Agreement be- 
tween the HI and stellar derived mass is roughly consis- 
tent when considering both the receding and approach- 
ing sides of the HI data. Differences most likely result 
from the uncertainty in assuming a velocity anisotropy 
for the stellar circular velocity reconstructions, and pos- 
sibly any non-circular motions not taken into account in 
the HI analysis. A weighted best fit isothermal sphere 
model based on the six stellar bins is shown in magenta 
and an NFW mass model in green. The rotationally de- 
rived mass at the half light radii (1656 ± 49 pc) from 

4.3 ± 0.3 x 10 8 M m , and 



the stellar data is Mi 



half 



Mi as t — 8.5 x 10 M at the last measured stellar bin 
(2865 pc). Assuming the best fitting enclosed mass pro- 
file from the HI study of |Kepley et al.| ( |2007[ ), one finds 
masses of M ha if ~ 3 x 10 s M Q , and M iast ~ 7 x 10 8 
Mq from the gas rotation at those same radii - in close 
agreement with the values from the stellar rotation. The 
NFW and isothermal sphere models that are fit to the 
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stellar data imply halo masses (taken to be the mass in- 
terior to a radius where the enclosed density equals 200 
times the critical density) of 8.9 x 10 9 M Q and 2.6 x 10 10 
M respectively. We show these values plotted versus 



the known stel lar mass of WLM (1.1 x 10 7 M Q ; [Jack- 
son et aL||2007 ) in the inset of Fi gure 15 along with the 
M* — Mhaio relation of Guo et al. (2010). The isothermal 
sphere model and NFW model are offset from each o ther, 
but both ro ughly consistent with the relation from Guo 
et al. (20101. However the sensitivity of the derived halo 
mass to the rotation signature at small radii in WLM 
prevents us from discussing the implications of this in 
any further detail in this work. With the above v alues 
and assum ing a V band luminosity of 5.2 x 10 7 L (Ma- 
teo||1998 ' ), we find a mass to light ratio of M/L = 8 and 
29 Mq Lq 1 at the half light radius, and 4.5 kpc (roughly 
the extent of the stellar and gaseous body of WLM) re- 
spectively. 

5.2. Dispersion Based Mass 

An alternative estimate of the dynamical mass can be 
taken from simple analytic expressions relating the mass 
to the line of sight velocity dispers ion at a characteristi c 
radius, as in the recent works of Walker et al.| ( |2009[ ), 
Wolf et al.| ( |2010[ ), and |Amorisco fc Evans] ( |201lfl inde- 
pendently, the studies identified a characteristic radius 
where velocity anisotropy effects are minimized. With 
these formulae the enclosed mass is co mputed in accor- 
dance with simple virial estimates. The Wolf et al. ( 2010 \ 
relation requires a luminosity weighted velocity disper- 
sion, which may introduce bias when applied to compos- 
ite populations wh ich have spatially differen t luminosity 



distributions (c.f., Amorisc o fc Evans 



2011). Therefore 



we use the relation of Walker et al. (2009) who found 
M(rh) — nr h a 2 , where \i = 580M Q pc km~ 2 s 2 . Being 
easily computed this is used frequently to derive masses 
in the dSphs and UFDs of the Local Group, however 
there are several points of caution to be noted in this 
formula's application to WLM. Any of these dispersion 
based mass estimators assume a system that is spher- 
ically symmetric, without rotation, and with a flat ra- 
dial velocity dispersion profile - all conditions which are 
violated to some degree in WLM. Nevertheless it is in- 
structive for the time being to compare the dispersion 
based mass to other estimators. We leave discussion of 
the biases that violation of these assumptions creates to 
our forthcoming paper on the dynamical modelling of 
WLM. With these caveats we proceed and compute the 
average velocity dispersion within using the rota tion 
subtracted velocity dispersion pr ofile shown in Figure [TTj 
With this the half mass from the Walker et al. ( 2009 ) re- 
lation is M ha if = 2.1 ± 0.3 x 10 s M . In Table 3 we 
summarize the mass estimates for WLM for ease of com- 
parison. 

With the enclosed mass at the half light radius it is 
now instructive to compare WLM to the half mass es- 
timates of Local Group dSphs. In Figure 16 we show 
that WLM lies slightly on the mo re massive side of th e 
sequence of M hal f vs. shown in Walker et al. ( 2009 ) . 
The other dSphs of the Local Group that lie off the rela- 
tion in this direction are the isolated Cetus and Tucana 



radius for a galaxy that has lost various percentages of its 
dark and luminous mass due to tidal stripping. This may 
provide an order of magnitude estimate for the impact 
of stripping on WLM, but the models assume a dwarf 
galaxy with a King profile and minimal rotation, so only 
a qualitative comparison is possible. While WLM's loca- 
tion on this diagram is not evidence itself for an isolated 
history, we can say that if it were to go through signifi- 
cant tidal evolution it would likely transition down into 
the sequence of the diagram occupied by the majority of 
the dSphs. 

A more appropriate comparison may be to consider 
WLM with respect to the dwarf elliptical (dE) satellites 
of M31. Studies of NGC 147, NGC 185, and NGC 205 
have shown them to be of similar total mass and metal- 



licity as WLM ( 


Bender et al.|1991 


De Rijcke et al.|2006 


Geha et al.||201 


J) although with hi 


gher luminosities anc 



( |Mateo|1998) . Analyses of the v elocity fields of NGC 147 

X„A \\/ V < 1 Jc l, — „„j;; l,,, \r*„u„ „-t „i I /lorn nh „,.„ 



and JNIGC 185 at large radii by Geha et al. (2010) sug- 
gest that these objects have similar kinematics as WLM, 
with all three exhibiting well defined rotation. However 
given the larger axial ratios of NGC 147 and NG C 185, 



their (V/a) * values are 0.95 and 0.91 respectively (|Geha 
eTalTpOlOt compared to (V/a)* ~ 1.1 for WLM 
tore while the dEs are more similar dynamically to WLM 
than the dSphs, the dEs likely have stronger anisotropy 
amoung their stars compared to WLM based on this ra- 
tio. Additionally the mass to light ratio (M/L) of WLM 
is larger than the dEs despite having similar total masses 
at comparable radii. Specifically, while the M/L ratios 
for the thr ee dEs are all consiste nt with M/L = 4 in 



solar units ( De Rijcke et al~| 2006), WLM's ranges from 
16 < M/L < 30 over the half light radius to the last 
me asur ed point. From our circular velocity curve in Fig- 
ure [15] we estimate WLM's dynamical mass at two ef- 
fective radii to be M(2R e ) ~ 1 x 1O 9 M . This is very 
near to the mass of NGC 147 and NGC 185 at simi- 
lar radius (5.6 x 10 8 ,7.2 x 10 8 ; |Geha et al.[|2010|), bu t 
slightly lower than NGC 205 (2 x 10 a ; |Geha et al.|2006 ). 
While WLM dynamically resembles ISGC 147 in partic- 
ular, it would seem that the relative dark matter content 
of WLM is somewhat higher than all of the Local Group 
dEs. However it is interesting to note that a conversion of 
WLM's remaining gas mass to stellar mass would result 
in a (M/L) ~ 4, consistent with the dEs. 

We note that the derived mass for WLM is above the 
possible threshold for gas retention found in the dwarf 
galaxy simulations of Sawala et al. (2011). That paper, 
along with rece nt simulation and observa tion work in the 



last two yea rs (Governato et al. 2010 



et al. 2010), indicate that the globa 



dSphs. Tidal evolution tracks of Penarrubia et al. ( 2008 1 
are overlaid, showing the expected evolution in mass and 



Sanchez- Janssen 
structural and 
dynamical properties of dwarf galaxies with gas are 
more strongly impacted by the initial halo mass, and 
their subsequent evolution driven by internal feedback 
or secular evolution, rather than tidal effects from a 
larger host galaxy. However environment likely still 
plays some role given that WLM and the dEs (which 
both appear to have comparable total masses within two 
effective radii) have drastically different gas fractions 
and effective radii - presumably due to the fact that 
the dEs reside within the virial radius (< 200kpc) of 
M31 where ram pressure and tidal stripping is expected 
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to be efficient. Differences such as these make it even 
more important to characterize the evolution of isolated 
galaxies, and the analysis in this work suggests that 
WLM represents one possible outcome of dwarf galaxy 
evolution in the absence of strong environmental forces. 



6. DISCUSSION 

In the last sections we have presented data on the 
structure and dynamics of WLM, which shows a thick 
gaseous and stellar configuration, as well as equal 
pressure and rotation support in the stellar kinematics. 
Constraints based on the current position and velocity 
of WLM with respect to the Local Group barycentre 
indicate that WLM could only have had at most one, 
if any, pericentre passage(s) and spent most of its 
history in isolation, with minimal tidal interactions 
in its lifetime. Thus we might ask how the extended 
structure of the gas and stellar disks/spheroids, and 
relatively high p ressure support in the s tars came about? 
As di s cussed in Brook s~et al.| |2007[ ) : [Opp e nheimer fc | 
Dave (20081: Sanchez- J anssen et al. d2010k 



et al. 



12008); Sanchez- J anssen et al. (2010); Govcrnato 
( 2UTU l; |Roychowdhury et al.| ( |2010p ; |Sawala et al. 
3nTT| thi effect of stellar feedbicTr~(SFR, HNe) is 



increasingly important in galaxies as their total mass 
decreases. Other internal processes that can work 
to transform the orbits and structure of stellar and 
gaseous components include global disk inst abilities, 
bars, and giant molecular cloud (GMC) heatin g (Spitzer 
fc Sch warzschild |1951[ |1953| |Sotnikova & Rodionov 



2003p . These effects may also increase with time in a 



system with ongoing SF like WLM. We can gain insight 
into these processes by studying the change in velocity 
dispersion over the course of WLM's history. 



6.1. Velocity Dispersion Evolution 

In Figure K7\ we show the velocity dispersion in stellar 
populations of different ages. There is a trend of increas- 
ing velocity dispersion with increasing stellar age, and at 
all ages the stellar dispersion is above that of the HI. It is 
useful now to discuss how exactly the stellar dispersions 
may have increased ~ 10 km s _1 higher than the neutral 
gas velocity dispersion. One interpretation of Figure 
is that the stars of all ages are heated at some rate, an< 
originated with an initial, constant dispersion - which in 
our c ase agrees with the m easured HI velocity dispersion 
from 



Kepley et"aL] ( |2007| ) (dotted line^ Alternatively, 



stars of different ages could have formed with different 
dispersions depending on the ISM kinematics at the time 
of their birth. Given that it is possible both interpreta- 
tions are working in combination, we now briefly discuss 
a few processes that might be expected to contribute to 
heating of the stellar populations in WLM. 

6.1.1. Tidal Heating 



Mayer et al.| ( |2001b I and |Mayer et al.| (|2006j) showed 
that tidal effects of a dwarf galaxy passing through peri- 
centre on its orbit about the Milky Way may trigger bar 

2 The dispersion of the gas that stars formed out of may be 
even colder than the HI dispersion if one considers the formation 
occurring in the cold molecular phase 



formation, which later undergoes a bending instability 
and then works to increase dispersion in the stellar com- 
ponents. However, this bar formation is not expected to 
occur (at least to a strong degre e) in populations that ar e 
already heated, or born thick ( |Kazantzidis et aL 2011). 
Recent simulations that use dlrr progenitors modelled 
with more realistic S F and feedback effe cts and inject 
them into the orbits of Mayer et al. ( 2006 1 seem to show 
that the structural and dynamical transformation is also 
possible through direct tidal interactions between the 
stars of the d warf galaxy and the potential of the MW 
( |Mayer|[2011[). This c o nclusi on appears to be at odds 
with IPeriarrubia et al. (2008), who showed that disper- 
sion decreases at all radii as tidal stripping progresses - 
a trend that goes in the opposite direction than what is 
seen in Figure [17} However, if new stars continue to form 
as the dwarf is tidal ly stripped of mass (a scenario not in- 
cluded in Penarrubia et al. ( 2008[ )), perhaps both models 
could lead to lower velocity dispersion at younge r ages . 



We also note that even the simulations of Mayer] (2011 1 
that produce hot spheroids through direct tidal heating 
only place their dwarfs on orbits that have apocentres of 
r apo — 250 kpc. The implied apocentre of r apo > 1 Mpc 
for WLM (Figure PL] and orbital time of > 11 Gyr, lim- 
its it to having haaat most one pericentre passage, and 
would seem to rule out such a tidal scenario for WLMs 
heating curve. 

6.1.2. ISM Pressure Floor 

Alternatively, Figure [17] could be interpreted as show- 
ing that the oldest stars formed from birth in an orig- 
inally hotter configuration. Simulations by |Kaufmann 
et al. (20071 found that in lower mass halos, an ISM 



temperature floor could produce initially thicker dwarf 
galaxies, as gas pressure support becomes increasingly 
dominant in lower mass halos. If this ISM temperature 
floor decreased with time, perhaps a fraction of the evolu- 
tion of the stellar velocity dispersion could be attributed 
to the stars forming out of a hotter ISM at earlier epochs. 
In the "cold flow" picture ( Keres et al.|2005||Brooks et al. 
2009 ) , this galaxy should at all times have accreted gas 
"cold" (i.e., on the free fall time), but at high redshift it 
is likely along filaments, while below z ~ 1 — 2 it would 
be more spherical accretion. In this case the filamentary 
accretion perhaps could drive turbulence in the ISM to a 
strong degree - which would naturally settle with time as 
the accretion became more spherical. Explicit evidence 
for accretion is difficult to obtain, but the metallicity dis- 
tribution function for the WLM stars does show a devia- 
tion from a closed box chemical evolution model, and in 
addition the spread in metallicity shows a slight increase 
as a function of age (Leaman et al., 2012b). 

6.1.3. Giant Molecular Cloud and Dark Matter Subhalo 

Heating 

The merger history of an object like WLM may also 
be important in dictating the evolution of its velocity 
disper sion and vertical thickness w ith time. Abadi et al. 
(2003) and |Robertson et aL ( |2006[ ) looked at the relative 
contribution of mergers of varying masses and gas frac- 
tions. The latter study noted the importance of gas rich 
merging as a way to maintain rotation su pport in galax- 
ies, while the work of Abadi et al. (2003) showed that a 
hot, thick stellar disk often forms early on, around the 
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time of the last major merger. However both studies 
could not include the high resol ution stellar feedback ef - 
fects shown to be important in Governato et al. (2010), 
and were simulating galaxies that were several orders of 
magnitude more massive than WLM. 

The case of heating by giant molecular clouds (GMCs) 
usually assumes that the clouds are occupying a rela- 
tively thin region of a galaxy and the encounter velocity 
is dictated by the peculiar velocity of th e stars in a 
rotating frame. Therefore cx ^ (c.f. 



2001 1, and the heating is not particularly efficient as the 



Fuchs et al. 



rate drops as the random velocity component increases. 
This is due to the reduced encounter probability, as the 
stars spend more time outside the GMC "layer" as the 
dispersion increases. The blue solid line in Figure [17] 
shows a canonical relation of the form: 



a{t) cx {t-t Y 



(4) 



The maximum heating that is expected analytically and 
from simulations o f giant molec ular cloud collisions alone 
is olqmc < 0.25 (Lacey 1984). While the mass spec- 
trum of GMCs in WLM is not precisely known, there is 
observational evidence to suggest t he properties are no t 
strongly dependent on host galaxy ( Bolatto et al.| 2008 ) . 
This maximum GMC contribution (indicated by the blue 
line) produces a curve that comes quite close to repro- 
ducing the observed stellar dispersion values. 

An additional source of heating may be present due 
to interactions between dark matter subhalos and the 
stellar component of WLM (c.f. 
These perturbations would have 



Moore et 
higher 



al. 1999) 



encounter 

velocities and expected heat ing would be more or less 



isotropic (Fuchs et al. 



120011). A simple calculation for 
the amount ot heating that a host galaxy of halo mass 
Mh undergoes from minor mergers with dark matter 
subhalos of mass M s can be derived assuming that the 
rate of heating for a subhalo to deposit energy into the 
disk component of the galaxy is: 



8E d = AE d 

at 



df 



where, t d f is the dynamical friction timescale. 



tdf - 



1.17 



'M h 



ln(M h /M s ) \M S 



r_h 



(5) 



(6) 



For a given interaction, the amount of energy deposited 
into the halo component relative to the disk comp onent 
of WLM scales as rj cx 2.3M h /M d ~ 27 for WLM pbth 
& Ostriker 1992). Thus the actual energy input into 
the disk will not be a particularly large fraction of the 
satellite kinetic energy. Using the expected unevolved 
subhalo mass function which is shown to be robust to 



halo size and over a range of redshifts (Giocoli et al. 
20081): 



dn 



51n(M s /M h ) 



= A 



( M a \-p r / 



fM h 



M s 
-fMh 



(7) 



Mo et al. (12010) show that the total heating rate by dark 
matter subhalos is: 



dE d 

at 



X 



n{M s \M h )M*dM s , 



(8) 



where the expected number of subhalos for a given pri- 
mary halo can be further simplified to n(M s \M} l ) oc 
M~ 18 . In order to equate this heating rate by the dark 
matter halos to an observable, we make the assumption 
that the kinetic energy deposited in the disk by a single 
event is: 

AE d = (l/2)M d a 2 z {R). (9) 

The total heating rate due to dark matter subhalos from 
masses m m i n to m max can then be rewritten as: 



da^ 2y / (3)M 1 - s 
dt ^ 1.2M d 



(10) 



Therefore the expected total heating rate by dark mat- 
ter substructure can be roughly estimated by fitting our 
velocity dispersion-age data with a function of two vari- 
ables; the maximum subhalo mass (m max ) and the WLM 
disk mass (M d ). We assume equipartition of energy in 
all velocity directions, which has been shown to be ap- 
propriate for satellites with circular orbits of various in- 
clinati ons, and an isotropic M axwellian distribution of 
orbits (Toth & Ostriker 1992). Similarly, assuming an 
average dispersion within lr e for an age bin is justified, 
as the scale affected by the energy inject ion of satellites 
has b een shown to be at least that large ( Sellwood et al. 
19981. Therefore it is a reasonable assumption for this 



simple model that disk stars at nearly all radii are heated. 
In red we overlay two examples for a heating model based 
solely on DM interactions. Both curves have a maximum 
perturber corresponding to m max — 5 x 10 8 M Q , or ~ 5% 
of the total halo mass of WLM. The dashed line uses a 
disk mass of 1 x 10 8 M , close to the present day value, 
while the solid line assumes a disk mass of 8 x 10 5 which 
is approximately the di sk mass of WLM 1 1 Gyr ago usin g 
the scaling relations of Dutton & van den Bosch (2009). 
Given the assumptions discussed in this paragraph, both 
curves are strong upper limits on the amount of heating 
substructure would predict. Despite this, the fact that 
the energy deposited in the disk scales as the mass ra- 
tio of WLM's disk to halo, prevents strong contribution 
from this effect. Even computing substructure models 
with an unphysical value of rj = 1 does not match the 
observed dispersion evolution. This is in good agree- 
ment with sem i-analytic and n-bod y simulations of disk 
heating such as |Benson et al. (2004), which found similar 
inefficiencies for substructure heating. 

In that study the authors showed that for lower mass 
primary galaxies, the contribution of heating from GMC 
encounters will dominate relative to substructure interac- 
tions. As an additional independent check, we compute 
th e heating/thickening predicted by Equations 11 and 12 



of Benson et al. (200 1 i. expressed as a ratio of vertical 



to radial scale lengths, h = z x /r d . With modest assump- 
tions for the parameters of those equations we find the 
same trend - GMC heating can produce a system with 
0.32 < h < 0.57, close to the observed inclination cor- 
rected axial ratio of WLM (see §4.4). Using the same 
range of disk and perturber masses in our Equation 10, 
we again find that the substructure mechanism is not 
enough to produce the observed vertical scale height - 
with only 0.02 < h < 0.09 found at the end of 12 Gyr 
in those models. However recent simulations have shown 
that disk heating by dark subhalos may be more effec- 
tive if higher mass ratio accretion events on plunging 
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radial orbits are considered, however the probability of 
such events occurring is still thought to be low (Tjitske 
Starkenburg, private communication). While there are 
obvious degeneracies that prevent us from accurately 
knowing the exact portion of the heating that came from 
GMC or DM substructure encounters (especially given 
the uncertain merger rate as a function of time), com- 
bined they are clearly of the right order of magnitude to 
explain the increased velocity dispersion and thickening 
in WLM. Most importantly these processes can occur in 
a galaxy in isolation, and don't require influence from 
the MW or M31 to modify a galaxy. 

Additional secular processes may work to heat the 
disk in the vertical direction, including long scale length 
global disk instabilities, as well as transient spiral o r 
overdensity features (c.f., Sotnikova & Rodionov 2003 1 . 
However, as those features may only last 1 Gyr, we can 
not comment on their impact directly. Although we have 
not ruled out a scenario in which gas dispersion decreases 
with ti me (§6.1.2), the dispersion evolution presented in 
Figure \T% when coupled with WLM's isolation, strongly 
suggests that the gaseous and stellar populations became 
heated and thickened through continual internal effects 
such as GMC and substructure interactions, rather than 
tidal transformations triggered by the MW. 



6.2. Comparing WLM to Evolutionary Models 

How do these observations fit in with the larger picture 
of dwarf galaxy formation and evolution? The structure 
and dynamical evolution of this isolated dwarf galaxy can 
be examined in perhaps one scenario that is consistent 
with th e findings presented i n this paper. 

From [Sawala et ah (2011) we see that the ability of 
a dwarf galaxy to retain gas until the present day may 
be due in most part to the total halo mass of the dwarf 



- with a threshold of approximately 1 x 10 9 M Q below 
which the galaxies are less likely to retain gas. Gas rich 
galaxies in that simulation also tended to have a late 
infall (z < 0.2) to the group, and pericentre distances 
that were greater than 300 kpc. Because of this, dlrrs 
could have gas over a larger percentage of their history, 
and experience ongoing SF and SNe - internal feedback 
that contributes to heating the dwarf, creating a slightly 
puffier syste m. Additionally, simulations by |Schroyen"| 
et al. (20111 show that dwarf galaxies with high angu- 
lar momentum tend to have SF continuously over the 
full radial extent of the galaxy - compared to non rotat- 
ing systems in which the gas can efficiently fall to the 
center and undergo extr eme centralized SF and blowou t 
events. As discussed in Robertson & Kravtsov (2008), 
if SF is regulated by the ability of the gas to settle into 
the molecular phase, a thick, hot ISM would naturally 
result in low efficiency SF, in line with the high gas frac- 
tion seen for WLM today. Therefore higher mass, rotat- 
ing dwarf galaxies would more likely be the type where 
we would expect to see signatures of extended internal 
feedback effects - namely a thick stellar configuration. 
These systems, if isolated, will be solely subject to inter- 
nal feedback effects such as SF, GMC interactions and 



secular glo bal disk instabilities (c.f. Sotnikova & Rodi- 
onov| 2b03 ) , which only produce mild transformations in 
their structure and dynamics. 



Thus we could have a scenario where the initial mass 
dictates th e probability of gas retention (i.e., 



et al. 



Sawala 



2011), which in combination with angular momen- 
di 



turn (Schrdyen et al. 20111 Koleva et al. 2009), sets the 



level of internal teedback driven heating. In this picture 
WLM could be an example of an isolated dwarf irregular 
at high enough mass and isolation to allow for ongoing 
SF driven feedback which maintains a dynamically hot 
stellar structure, low SF rate, and high gas content. 

7. SUMMARY AND CONCLUSIONS 

We have presented the results of a VLT and Keck II 
spectroscopic survey of 180 RGB stars in the isolated dlrr 
galaxy, WLM. Broadband optical and infrared photom- 
etry and HI radio observations allowed for the study of 
the structure and dynamics of stellar and gaseous popu- 
lations in a galaxy at the outskirts of our Local Group - 
a unique oppourtunity in resolved spectroscopic studies 
of dwarf galaxies. 

The main results from our analysis are: 

• Given the high inclination, the apparent axial 
ratios of both the gas and stars indicate that 
WLM has a thickened structure (q) = 0.39 - 0.57. 
However the data is not detailed enough to dis- 
tinguish between a thick disk, or flattened oblate 
spheroid. This thickened structure is consistent 
with ob servations of statistically large samples of 



dwarfs ((S anchez- Janssen et 
hury et al.||2UlUp 7 



al. 



2010 Roychowd- 



• The ellipticities and position angles of the HI, 
young C-stars and evolved RGB stars are nearly 
identical. The photometric and kinematic position 
angles of the gas and stars are also closely aligned 



gas 



2°). 



• We measure a maximum rotation in the stars 
of V rot = 17 km/s, which after correction for 
asymmetric drift approaches V ro t = 36 km/s - 
in excellent agreement with the HI. The velocity 
dispersion profile reaches a maximum of a ~ 17 
km/s and declines with radius, unlike the flat 
profiles seen in dSphs. The ratio of rotation to 
pressure support, V/a, is of order unity out to 
1.5r e , where it slowly rises to 2.5 - in contrast to 
the HI, which shows V/a ~ 7. 



The stellar velocities show V/ a ratios that decrease 
with age, and the velocity dispersion grows with 
age at a rate that can be explained through pri- 
marily GMC heating and substructure interactions 
- with the former dominating. The evolution of 
o~(t) is therefore consistent with internal feedback 
mechanisms, rather than external environmental 
effects such as tides. 



• Using rotationally derived mass estimates we com- 
pute a mass at the half light radius (r^ = 1656 ±49 
pc) of M hai/ ~4x 10 8 M . 
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• Put together, the observations paint a picture in 
which WLM has spent the majority of its evolution 
in isolation, with its velocity dispersion and struc- 
tural evolution governed primarily by conditions at 
formation and subsequent internal feedback effects, 
rather than tidal events. 

This departure from dwarf irregulars as thin rotation- 
ally supported cold disks may have implications for the 
timescales required for tidal transformation scenarios, 
and for analysis of the statistical properties of dwarf 
galaxy populations on the group scale. WLM represents 
one look at how only the baseline secular processes in a 
galaxy will modify it over ~ 12 Gyr. Further studies of 
dlrr galaxies of differing masses and isolation levels will 
continue to disentangle to what degree internal feedback 
and environmental evolution play in forming the distinct 
classes of Local Group dwarf galaxies we see today. 
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Fig. 1. — Distance from M31 and the Milky Way for Local Group dwarf irregulars (blue), dwarf spheroidals (red), and transition dwarfs 
(green). Shown are the projected galactocentric standard of rest (GSR) and Local Group standard of rest velocities for WLM, Eviden t 
is WLM's large isolation from the two massive spiral galaxies, as well as other dwarf galaxies. Coordinates are taken from [Mateo| |l998[ l . 
Shaded areas correspond to the approximate virial radii of the Milky Way and M31. 
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Fig. 2. — Portion of an I band image of WLM from the MOSAIC-II camera at the 4m Blanco Telescope at CTIO (Leaman et al., in 
prep.). The total image is approximately 36' X 36', equal to the coverage from our INT WFC photometry. This Figure shows a zoomed in 
region of 22' X 26', with North being up and East to the left. The relative locations of the VLT FORS2 (blue), and Keck II DEIMOS (red) 
spectroscopic fields are shown. Magenta ellipse marks the half l ight radius assuming an ellipticity of e = 0.55 and photometric position 
angle of 179°. The two A- type supergiants from l V enn et al. I d2003[l are located approximately in the centre of our lower FORS2 field, along 
with the HII regions from |Hodge fc Miller| ( |l995| l, and the B supergiants from [Bresolin et al.| l [2006| . 
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Fig. 3. — Example DEIMOS co-added spectra for three member stars in our sample of varying magnitudes, plotted in arbitrary units of 
flux and continuum normalized and binned by a factor of 4. Visible are the three strong Ca II triplet lines at A ~ 8498, 8542, and 8662A 
shown by the red dotted lines, as well as the DEIMOS chip gap (blue regions). The signal to noise per angstrom for the spectra are ~ 30, 
19, and 15 respectively. 
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Fig. 4. — Heliocentric velocity as a function of major axis position for the 180 member RGB stars in WLM, as well as the non-member 
stars. Extremely large absolute velocities are typically stars with very low signal to noise which rendered clean sky subtraction and radial 
velocity derivation difficult. The mean random uncertainty on velocity for the member stars is 4 km s _1 (smaller than the symbols in most 
cases). Due to the large distance off the plane of the galaxy (6 = —74°), few Milky Way stars contaminate our spectroscopic catalogue. 
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Fig. 5. — WLM spectroscopic sample in tangent plane coordinates on the sky. HI velocity data from |Kepley et al.| l |20(j7] l is shown as 
ontours, and the stellar velocities of RGB members are colour coded to the same velocity scale. Dashed line shows the photometric minor 



contours, 

axis, and the dotted line the kinematic angle of bisection determined as discussed in text. 
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Fig . 6.— INT WFC V,I photometry, with strict RGB cut, presented with JHK UKIRT/WFC AM photometry of WLM ( |Tatton et aT] 
|2010 |. The metallicity d ependent cut between C-rich AGB stars (blue box) and the O-rich AGB stars for WLM (light blue box ) is taken 
lfom |Tatton et a.l.| J2010| | . 
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Fig. 7. — Geometrical radial surface density profiles for the full V,I (black), strict RGB (red), and C-star (blue) stellar populations. 
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Fig. 8. — Density contours for the stellar and gaseous populations of WLM in the tangent plane. Contour levels start at 2a above the 
background and each subsequent contour level is n times higher than the last, where n = 1.2. 
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Fig. 9. — Derived ellipticity and photometric position angle as a function of major axis radius for the various populations. Black is the 
full INT catalogue of stars, r ed the strict RGB pop ulation, blue the C-rich AGB stars, and mage nta the HI gas. Dashed l ines show the 
measured disk ellipticity from | Abies & Abies ( 1977|l in black, and extended RGB halo ellipticity of |Minniti fc Zijlstra ( 1997|l (red). 
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Fig. 10. — WLM stellar and HI rotational velocities (V^i — V aya ) as a function of major axis offset. Small black dots represent the 
individual stellar data points, and the large circles the binned values with error on the mean shown. Results of the asymmetric drift 
correction to the stellar rotation is shown in blue. Also noted are the relative rotation over pressure support for the raw stellar and gaseous 
kinematics. 
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Fig. 11. — Raw stellar rotational velocity (top), along with the optimal derotation achieved through a single fit through the major axis 
velocity profile of WLM. Third panel shows velocity dispersion profile before (blue circles) and after rotation is subtracted - with the 
artificial enhancement of <j v at large radii removed in the latter case. Bottom panel shows ratio of rotational to pressure support, with 
dashed line showing the expected value for an oblate galaxy flattened by isotropic rotation to the ellipticity of WLM. 
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Fig. 12. — Kinematic position angle diagnostic plots for the HI (top) and stellar data (bottom). The velocity difference (sec text) between 
the stellar sample on either half of a bisecting axis is plotted as a function of the bisection angle. Here the bisection angle of zero degrees 
corresponds to the minor axis of WLM. Overlaid is a function of the form aotanh(ai(0 — a2)). 
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Fig. 13. — Range of expected intrinsic thickn ess for an observed p rojected axial ratio. Shown are curves for the mean, upper and lower 
inclinations from the HI tilted ring analysis oflKepley et al. M20071. Horizontal dotted lines represent the range of axial ratios from the 
isopleth analysis presented in Figure [9] with the horizontal dashed lines the mean values for the stellar populations. Vertical dotted lines 
indicate the corresponding range of go values given the inclination limits from the HI analy sis. Vertical dashdot line indicates the intrinsic 
axial ratio assuming the measured axial ratio from the minor and major axis fits of Figure [14] 
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Fig. 14. — Comparison of m ajor axis radial scale lengths versus vertical thick disk sc ale heights. Black dots are data from the work of 
|Yoachim fc Dalcanton (2006 ), and M31 and MW data taken from [Collins et al.| l |20lT| | and references t herein. Inset show s vertical scale 
heights m parsecs as a function of V c i rc for the same data pl us 6 resolved stellar pop ulation studies from |Seth et al.| (|2005|> ■ Black line is a 



fit of the form Z 0:thic k = (1400pc)(- 



— ) adopted by 



' lOOkms 

shows the effect of the inclination correction as discussed in the text 
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I 2006 ) for their galaxy sample. Open circle for WLM 
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Fig. 15. — Enclosed Mass as a function of geometrical radius. Black squares are rotationally derived mass estimates using the asymmetric 
dr ift corrected stellar v elocities at various radii. Black lines show the HI rotation curve isothermal sphere estimates from the best fit models 
oflKepley et al.l (120071 to the receding and approaching sides of the galaxy. Cyan curve is the best fitting isothermal sphere model fit to 
the six stellar data points (see text). Green curve is NFW mass model fit to the same points. Inset shows stellar mass as a function of the 
implied halo mass from both the SIS and NFW models, with comparison to the relation (black line) from |GucT ct al. ( 2010j and references 
therein. 
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Fig. 16. — Half mass as a function of half light radius for Local Group dSphs, and WLM (blue). Overlaid are the expected tidal evolution 
tracks of Pcnarrubia et al. 2008 for various total stripping values of the initial luminous mass. WLM lies on the isolated side of the relation 
which is consistent with its implied separation and orbital history (see §1.1). 
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Fig. 17. — Velocity dispersion as a function of age in the stellar components of WLM. Random error on the mean and systematic errors 
for each age bin ar e shown as the solid and dash-dot error bars respectively. (Youngest data point (t = 0.01) is calculated from the 
supergiant study oflBresolin et al.M2006t. Blue dashed line indicates the maximum expected heating from simulations of giant molecular 
cloud heating, and red solid and dasned line represent an upper and lower contribution from DM substructure heating (see text). Composite 
curves including both effects added in quadrature is shown in magenta. 
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TABLE 1 
WLM Properties 



Quantity 



Value 



Description 



Reference 



(a, 5) 

(1,6) 

E(B-V) 

Dmw 

e 

90 

zo 

vd 

Th 

h MA 

n 

yHI 

V hel 

M v 
M HI 

Mg 

M b 

Mb 
M tot 

Mtot 
L V 

r)9a 



T — — 

2irGI r HB 



(00 01 58, -15 27 45) 
(75.85, -73.63) 
0.035 (mag.) 
985 ± 33 kpc 

0.45 - 0.6 
~ 0.38 - 0.57 
705 ± 28 pc 
987 ± 30 pc 
1656 ± 49 pc 
1209 ± 23 pc 

178 ± 1° 

179 ± 2° 
-130 km s _1 

~ 6 
~ 1 
-14.1 (mag.) 
(6.3 ±0.3) x 10 7 M Q 
1.1 x 10 7 M 
0.86 
0.05 - 0.10 
17-30 M Q Lq 1 

232 ± 31 M Lg 1 
-1.45 ±0.2 dex 
-1.28 ±0.02 dex 



J2000 

Line of Sight Reddening 
Galactocentric Distance 
Eccentricity 
Intrinsic Axial Ratio 
Vertical Scale Length 
Radial Scale Length 
2D Projected Half Light Radius 
Major Axis Disk Scale Length 
Kinematic Position Angle 
Photometric Position Angle 
Heliocentric Velocity 
HI Rotational to Pressure Support 
Stellar Rotational to Pressure Support 
Integrated V Band Absolute Mag. 
HI Mass 
Stellar Mass 
Gas Fraction (M s = A/ H/ /0.733) 
Baryon Fraction (M b = M* + M g ) 
Direct Mass to Light Ratio (at r^ a if to r; ast ) 

Theoretical Mass to Light Ratio 
Photometric RGB Metallicity 
Spectroscopic RGB Metallicity 



Mateo 1 1998 1 



Gallouet et al. 



p975) 



|McConnachie et al,| ( 2005) 

Average 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 



Jackson et al. 


(2001 


Kepley et al. 


20071 



This work 
van den Bergh|(|1994} 



Kepley et al. 1 2007 \ 



Jackson et al. (2007) 



This work 
This work 
This work 

This work 



Minniti & Zijlstra 



(1997) 



This work 
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TABLE 2 
Repeated Measurements 



Matched Pair ID 



1 26940 

1 26939 

2 27429 

2 27428 

3 27926 

3 27925 

4 28220 

4 28219 

5 28395 

5 28394 

6 28986 
6 28985 



Vhel 


[Fc/Hl 

L / J 


-98.000 


-1.66 


-105.110 




-123.000 


-0.79 


-118.350 


-1.03 


-123.000 


-1.11 


-135.010 


-1.37 


-145.000 


-0.97 


-1333.610 




-139.000 


-1.24 


-124.940 




-112.000 


-0.76 


-104.620 





Note. — First entry of each pair is the 
FORS2 observation, second the DEIMOS. 



TABLE 3 

Mass Estimates 



Description 

Mass within 

Mass within from V* 

Mass within r* t from V* 

Mass within r^ t from approaching side HI SIS Model 
Mass within rfz I t from receding side HI SIS Model 
Mass within rf 11 . from stellar SIS Model 

last 

Mass within rf 11 , from stellar NFW Model 

last 

Mass within r 2 oo from stellar SIS model 
Mass within r'200 from stellar NFW model 



Value 

2.1 ± 0.3 X 10 s M 

4.3 ±0.3 x 10 s Mq 

8.5 ±0.3 x 10 s Mq 

9.4 ±0.5 x 10 s M 

2.2 ±0.1 x 10 9 M 
1.2 ±0.1 x 10 9 M 

1.6 ±0.1 x 10 9 M 
2.6 ±0.2 x 10 10 M Q 
8.9 ±0.8 x 10 9 M 



Notes 



Eqn. 1 of 



Walker et al 



(2009) 



Last measured stellar bin = 2.9 kpc 
Last measured HI point = 4.5 kpc 
Last measured HI point = 4.5 kpc 
Last measured HI point = 4.5 kpc 
Last measured HI point = 4.5 kpc 
r 200 (SIS) = 60.5 kpc 
r 2 oo(NFW) = 42.0 kpc 



